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The adaptive immune system relies on different cell types to provide fast and coordinated 
responses, characterized by recognition of pathogenic challenge, extensive cellular prolifer- 
ation and differentiation, as well as death. T cells are a subset of the adaptive immune 
cellular pool that recognize immunogenic peptides expressed on the surface of antigen- 
presenting cells by means of specialized receptors on their membrane. T cell receptor binding 
to ligand determines T cell responses at different times and locations during the life of a 
T cell. Current experimental evidence provides support to the following: (i) sufficiently 
long receptor -ligand engagements are required to initiate the T cell signalling cascade that 
results in productive signal transduction and (ii) counting devices are at work in T cells to 
allow signal accumulation, decoding and translation into biological responses. In the light 
of these results, we explore, with mathematical models, the timescales associated with T 
cell responses. We consider two different criteria: a stochastic one (the mean time it takes 
to have had TV receptor -ligand complexes bound for at least a dwell time, r, each) and one 
based on equilibrium (the time to reach a threshold number TV of receptor -ligand complexes). 
We have applied mathematical models to previous experiments in the context of thymic nega- 
tive selection and to recent two-dimensional experiments. Our results indicate that the 
stochastic criterion provides support to the thymic affinity threshold hypothesis, whereas 
the equilibrium one does not, and agrees with the ligand hierarchy experimentally established 
for thymic negative selection. 
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1. SUMMARY 

The binding properties of T cell receptors for self-pMHC 
(peptide -major histocompatibility complex) ligands 
are the basis for the selection in the thymus of a 
useful (MHC-restricted) and safe (self- tolerant) T cell 
repertoire. There exists a wealth of experimental sup- 
port for the following: (i) T cell receptors must be 
bound to their ligands for a sufficiently long time to 
initiate the T cell signalling cascade and (ii) T cells 
require T cell receptor signal accumulation, which will 
be translated into appropriate biological responses. We 
have made use of mathematical models to test two 
different hypotheses: (a) the timescale of a T cell 
response correlates with the time it takes to have had 
TV receptor -ligand complexes bound for at least a 
threshold dwell time, t, each and (b) the timescale of 
a T cell response correlates with the time a threshold 
number, TV, of TCRs must be occupied at equilibrium. 
We demonstrate that scenario (a) provides, for a given 
T cell receptor, a ligand hierarchy that agrees with 
that experimentally established for thymic negative 
selection, and an intuitive way to understand self- 
non-self discrimination of pMHC ligand. Our results 
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suggest that a very small number (fewer than 10) of 
cognate ligand molecules is sufficient to elicit a T 
cell response, which is consistent with the serial 
engagement model. 

2. INTRODUCTION 

The adaptive immune system relies on different cell 
types to provide fast and coordinated responses, charac- 
terized by recognition of the pathogenic challenge, 
extensive cellular proliferation (division) and differen- 
tiation, as well as cellular death. T cells are a subset 
of the adaptive immune cellular pool that recognize 
immunogenic peptides (non-covalently bound to MHC 
class I and class II molecules expressed on the surface 
of antigen-presenting cells (APCs)) by means of special- 
ized receptor molecules on their membrane. A (human 
or mouse) T cell expresses about 30 000 copies [1] of a 
T cell receptor molecule (TCR), whose ligands (usually 
referred to as antigens, in this context) are complexes 
composed of a peptide bound to an MHC molecule 
(pMHC). T cell receptors are both degenerate (a 
given TCR can recognize different pMHC complexes) 
and specific (single or point mutations to a pMHC 
complex can significantly alter recognition); yet 
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TCR-pMHC interactions have low affinities [2-6]. In 
vivo, pMHC complexes are expressed on the sur- 
face of APCs; each (human or mouse) APC displays 
around 100 000 different pMHC complexes on its sur- 
face [7-11]. It is through interactions with pMHC 
ligands that T cells become activated and differentiate 
into effector T cells, which elicit immune responses. 
Thus, in order to study the conditions under which T 
cells initiate a response, one needs to understand the 
dynamics of TCR-pMHC binding. 

T cells are derived from precursor cells that migrate 
from the bone marrow to the thymus, where they 
rearrange their receptor genes and generate a unique 
(clonotypic) TCR. In the thymus, immature T cells 
(or thymocytes) are exposed to an antigenic micro- 
environment orchestrated by APCs of different types 
[12] that subject thymocytes to a 'double test' by dis- 
playing a wide range of pMHC complexes, with 
peptides derived from household proteins (self-peptides 
or self-pMHC complexes). Owing to the stochastic 
nature of the gene rearrangements [13,14], some TCRs 
will not be able to recognize a self-pMHC ligand 
(TCRs that are not functional). Other TCRs will recog- 
nize it too well, and could give rise to mature T cells 
with the potential to generate autoimmune responses. 
Thus, the need for a thymic double test to check the 
functionality of a thymocyte (positive selection) and 
its state of tolerance, so that it does not recognize 
self-pMHC complexes with high affinities (negative 
selection) [15]. Thymic selection allows only 2-5% of 
all thymocytes to survive and migrate to the peripheral 
sites of the immune system (lymph nodes, spleen, etc.) 
[16,17], where they continuously recirculate via the 
blood, surveying the antigenic environment displayed 
once again by APCs. 

TCR-pMHC binding events determine T cell 
responses (survival, proliferation, differentiation or 
death) at different times and locations during the life 
of a T cell [18]. For example, Naeher et al. [19] have 
made use of a photo-affinity labelling system (that 
allows quantitative analysis of pMHC monomer bind- 
ing to TCR) to show that MHC class I restricted 
TCRs exhibit an affinity threshold during negative 
selection. In the light of these results, it is natural to 
consider the question [20]: how does the number of 
TCR-pMHC-bound complexes relate to the outcome 
of negative selection? Current evidence suggests that 
both the duration and the number of TCR-pMHC 
bindings play a role [21-23]. Valitutti et al. [24,25] 
have experimentally shown that a few hundred 
pMHC molecules can serially bind thousands of 
TCRs. Finally, Sykulev et al. [26] and Davis and col- 
leagues [27,28] have provided experimental data 
suggesting that a few agonist pMHC ligands (5-10) 
are sufficient to elicit a T cell response. This body of 
work provides support for the following two hypoth- 
eses: (i) TCR-pMHC engagement needs to be 
sufficiently long to result in productive signal trans- 
duction [29] and (ii) T cells can integrate signals; 
that is, counting devices are at work in T cells to 
allow signal accumulation, decoding and translation 
into biological responses [25]. These ideas have been 
explored by different groups: Palmer and Naeher 



have provided a biophysical basis for their affinity 
threshold for negative selection hypothesis [19,20], 
Dushek et al. [30] have developed a mathematical 'pro- 
ductive hit rate model', Chakraborty and colleagues 
[31-33] have developed statistical models of how T 
cells convert analogue inputs into digital outputs 
and van den Berg & Rand [34] have introduced the 
idea of a mean triggering rate. The first set of authors 
introduces the concepts of dwell time of individual 
TCR-pMHC complexes and productive TCR inter- 
actions, and compute the number of TCR-pMHC 
interactions required as a function of the TCR- 
pMHC complex half-life, for a given choice of dwell 
time and number of productive TCR interactions 
(see fig. 3 of Palmer & Naeher [20]). Current exper- 
imental evidence supports values of dwell time, r, of 
around 4 s and number of productive TCR interactions, 
TV, below 100 [20,35]. In this study, we make use of these 
ideas to provide a stochastic T cell response criterion 
based on a mathematical model of TCR-pMHC mole- 
cular interactions. The dynamics of a small number of 
TCR-pMHC binding events, as suggested by the 
experimental evidence mentioned earlier, is naturally 
described as a stochastic process, without the need to 
assume that TCR-pMHC association/dissociation 
kinetics has reached thermal equilibrium [36-41]. 

TCR-pMHC binding, and receptor -ligand binding 
in general, is described by reaction kinetics, assuming 
that the ligand is in solution and that receptors are 
membrane-bound on T cells [42,43]. The kinetics is gov- 
erned, for a given choice of receptor and ligand, by two 
rates, k + and that give the probability per receptor 
and per unit of time of a binding and an unbinding 
event, respectively [36-39,44]. The study of reaction 
kinetics for the system A + B ^ C is not only limited 
to the case of receptor -ligand interactions, but is of 
wide interest and has been applied to other problems, 
such as crystal growth, gene clustering, cellular metab- 
olism and catalytic efficiency of enzyme reactions 
[45-51]. From a thermodynamic perspective, it is natu- 
ral to assume that, if one waits long enough, forward 
(association) and backward (dissociation) reactions 
reach a steady-state or equilibrium [34]. Thus, in §3, 
we investigate an equilibrium dynamics model of 
TCR-pMHC association/dissociation. One candidate 
T cell response criterion that will be explored is that 
the timescale of a T cell response correlates with the 
time a threshold number, TV, of TCRs must be occupied 
at equilibrium, T N . However, given the experimental 
support behind the hypothesis that a few agonist 
pMHC ligands can suffice to trigger T cell responses 
[3,26,52] and Palmer's affinity threshold for negative 
selection [19,20,30], a stochastic approach seems more 
appropriate [38,39,53]. Furthermore and as discussed 
earlier, (i) sufficiently long TCR-pMHC engagements 
are required to initiate the signalling cascade, resulting 
in productive signal transduction [35], and (ii) T cells 
can integrate signals; that is, counting devices are at 
work in T cells to allow signal accumulation, decoding 
and translation into biological responses [25,23]. 

With this experimental and theoretical evidence 
in mind, we explore a different criterion, namely that 
T cell responses take place once a given number of 
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Table 1. A summary of binding data [19,20]. All constants have been introduced and denned in §5. 



cell type 


ligand 


K d (M) 




h/2 (s) 


k on (s^M" 1 ) 


K (s" 1 ) 




k- (s" 1 ) 


SP thymocyte 


4P at 37°C 


1.1 x 10" 


■7 


41 


153 691 


5.1230 x 10" 


10 


0.0169 


SP thymocyte 


4A at 37°C 


5.5 x 10" 


6 


0.8 


157 533 


5.2511 x 10" 


10 


0.8664 


SP thymocyte 


4N at 37°C 


5.8 x 10" 


■5 


0.08 


149 385 


4.9795 x 10" 


10 


8.6643 


DP thymocyte 


4P at 37°C 


8.8 x 10" 


■8 


39 


201 966 


6.7322 x 10" 


10 


0.01778 


DP thymocyte 


4A at 37°C 


2.2 x 10" 


6 


0.79 


398 161 


1.3272 x 10" 


9 


0.8760 


DP thymocyte 


4N at 37°C 


2.9 x 10" 


- 5 


0.23 


105 719 


3.5240 x 10" 


10 


3.0658 



TCRs (and not necessarily in a simultaneous way), TV, 
have been engaged with ligand for at least a dwell 
time, r, each. The stochastic criterion requires counting 
the number of productive bindings (a binding is pro- 
ductive if it lasts longer than the threshold dwell 
time, r, to capture the essence of the signalling cascade 
[39]). The first time at which this stochastic criterion is 
satisfied (a first passage time (FPT), as considered from 
a stochastic process point of view [54]) will be referred 
here, and in a biological context, to as the (first) time 
to signal initiation (TSI). We will derive expressions 
for its mean value, or mean time to signal initiation 
(MTSI), T(N, r), as a function of TV, the number of pro- 
ductive TCR-pMHC engagements and r, the dwell 
time, and its variance. We study the MTSI for different 
pMHC ligands, thymocytes and temperatures (both 
association and dissociation rates are temperature- 
dependent). We make use of recent data [19,55,56] to 
compare the equilibrium criterion versus the MTSI cri- 
terion, to explore the affinity threshold hypothesis and 
to confront two-dimensional and three-dimensional 
binding data with the stochastic criterion. 

The study has the following structure: §3 describes 
the main results and in §4 we explore the immunological 
consequences of the results. Finally, in §5, we provide 
the mathematical details of the stochastic model devel- 
oped, and how to derive the deterministic model as a 
limit of the stochastic model. We also provide analytical 
expressions for the mean and the variance of the TSI, as 
well as for the time to reach a threshold number N of 
receptor -ligand complexes, T N . 

3. RESULTS 

3.1. Receptor -ligand binding dynamics 

Our model of receptor -ligand binding is motivated by 
the experiments of Palmer and colleagues [19,20], 
measuring binding of soluble, monomeric pMHC 
ligands to live thymocytes from Tl TCR transgenic 
mice. The binding and unbinding reactions can be 
represented as follows: 




where the circle represents a free ligand pMHC and the 
open box an unbound TCR. 

We consider two different subsets of TCR transgenic 
Tl T cells (monoclonal TCR): pre-selection double 
positive thymocytes (DPs) and mature single positive 



thymocytes (SPs) [19]. DPs express on average 7V R = 
3000 TCRs and SPs express on average N R = 30 000 
TCRs on their surface. In the experiments, a panel of 
pMHC ligand complexes is used [19]. Here, we restrict 
ourselves to three, denoted 4P, 4A and 4N, whose bind- 
ing parameters are given in table 1. The negatively 
selecting ligand, 4P, has the highest complex half-life, 
ti/z, and lowest equilibrium dissociation constant, 
K d [44]. We note that the parameters K dl ty 2 and k on 
have been introduced and defined in §5.2. The ligand 
denoted 4N is positively selecting, with the lowest 
ti/2 and highest K d . The ligand denoted 4 A is called a 
threshold ligand [19]; it can act as a positively selecting 
or negatively selecting ligand depending on its con- 
centration [19]. For each ligand type and for a given 
temperature, the mathematical models require the 
association and dissociation rates for the TCR- 
pMHC interaction, k ± (see §5 and Lauffenburger & 
Linderman [44]). 

Our criterion is that T cell responses are initiated 
by discrete (stochastic) events and not by attaining 
the state of thermal equilibrium or steady-state. 
In order to support our criterion, in §3.2, we explore a 
stochastic model of receptor -ligand binding, and in 
§3.3 we analyse the deterministic limit of the 
mathematical model. The distinction between deter- 
ministic and stochastic approaches is not solely a 
mathematical one, but a choice that has its roots in 
the biochemical distinction between equilibrium and 
non-equilibrium dynamics. 

3.2. Stochastic criterion for T cell responses 

We first study the possibility that immunological 
responses of T cells are determined by the number of 
TCRs engaged for a minimum threshold time [20]. In 
order to explore this scenario, we develop a stochastic 
model of TCR binding to pMHC ligand, described in 
§5.1. We are interested in calculating the time it takes 
a T cell to reach (for the first time and not necessarily 
in a simultaneous fashion) N TCR engagements with 
pMHC ligands, such that each engagement lasts at 
least a dwell time r. This time defines the stochastic cri- 
terion of T cell responses and is denoted in what follows 
the first time to signal initiation (or TSI). The first time 
to signal initiation (or time to signal initiation, for 
short) is analogous to a FPT defined in stochastic dyna- 
mical systems [49,57]. We note the following implicit 
assumptions of the stochastic criterion: binding 
events, as well as unbinding events are considered 
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initial ligand concentration p(M) initial ligand concentration p(M) initial ligand concentration p(M) 



Figure 1. Mean time to signal initiation (MTSI) for a T cell according to the stochastic criterion. MTSI, T(TV,r), for single posi- 
tive thymocytes (SP), Tl TCRs at 37 degrees. Different panels stand for different pairs of values (TV,t): (a) for (10,1), (b) for 
(10,4), (c) for (10,8), (d) for (100,1), (e) for (100,4) and (/) for (100,8). Time units are seconds. The physiologically relevant 
range of initial ligand concentration is shaded grey, and the dotted lines correspond to a time scale of 1 min, 1 h and 
1 day, respectively. 



independent and identically distributed random vari- 
ables [58,59], and the times to bind, as well as the 
times to unbind are considered exponentially distribu- 
ted random processes (see §5.1 for mathematical 
definitions and choice of notation). In §5.1, we provide 
analytical expressions for the mean and the variance 
of the time to signal initiation. 

We now make use of the Tl TCR data described in 
Naeher et al. [19] for three different pMHC ligands 
(4P negatively selecting ligand, 4 A threshold ligand, 
4N positively selecting ligand) and summarized in 
table 1. In figure 1, we plot T(N,t) for SP thymocytes, 
as a function of the initial ligand concentration, for 
different values of TV and r (see §5.1 for mathematical 
definitions and choice of notation). 

From figure 1 (supported by equation (3.1)), we note 
the following properties of the MTSI. 

— For any value of the initial ligand concentration, and 
for any choice of TV and r, the stochastic criterion 
yields the shortest time to respond to the ligand 4P 
and the largest time to respond to the ligand 4N, in 
agreement with the experimental data of Naeher 
et al. [19]. Furthermore, 4A (or threshold ligand) 
displays in all cases (figure 1) an intermediate behav- 
iour. Thus, the MTSI provides qualitative support to 
the affinity threshold hypothesis introduced in 
Naeher et al [19]. Given an initial ligand concen- 
tration, a choice for (TV, r) and a time T for a T cell 
response, the positively selecting ligand will not be 
able to initiate a response within a immunologically 
relevant time. 

— There is no reversion of the hierarchy of ligands as a 
function of the initial ligand concentration. 



— An increase in either TV or r increases the value of 
MTSI. 

— As the temperature increases, the MTSI also 
increases (data not shown). 

We now proceed to present our results in greater detail. 

3.2.1. For single positive thymocytes the mean time to 
signal initiation is shorter when binding the negatively 
selecting ligand 

We have considered different scenarios, firstly for SP 
thymocytes [20,35]: (a) TV = 10 bindings and r = 1 s, 
(b) TV= 10 bindings and r=4s, (c) TV = 10 bindings 
and r= 8 s, (d) TV = 100 bindings and r = 1 s, (e) TV = 
100 bindings and r=4s and (/) TV = 100 bindings 
and r = 8 s. We compare the mean T cell response 
times (MTSI) of a positively selecting ligand (4N), a 
threshold ligand (4A) and a negative selecting ligand 
(4P), for varying initial ligand concentrations for the 
TCR system Tl. As seen in figure 1, the negative-select- 
ing ligand is characterized by the shortest MTSI in all 
cases considered. This is a result of both its higher k on 
and lower A; off , which means bindings occur more fre- 
quently and are more likely to last longer. This also 
accounts for the more pronounced differences in MTSI 
when either the threshold time, r, is increased from 1 
to 8 s or when the number of bindings required, TV, is 
increased from 10 to 100. 

A consequence of the stochastic criterion is that, the 
larger the number of bindings required, TV, the less impor- 
tant stochastic effects become. In this case, the coefficient 
of variation tends to zero as TV increases (see equation 
(5.17)). Immunologically, this is relevant as the criterion 
relegates TV to a role secondary to that of the dwell time r. 



J. R. Soc. Interface (2012) 



2860 A stochastic T cell response criterion J. Currie et al. 




10 -10 1Q -9 1Q -8 1Q -7 1Q -6 1Q -5 1Q -4 1Q -10 1Q -9 1Q -8 1Q -7 1Q -6 1Q -5 1Q -4 

initial ligand concentration p(M) initial ligand concentration p(M) 



Figure 2. Fraction of unbound T cell receptor at equilibrium, / eq , computed making use of the equations derived in §5.2. 
(a) Double positive thymocytes, (b) Single positive thymocytes. Squares are the experimental measured values provided in 
Naeher et al [19]. 




l o -5 H 1 1 1 1 1 1 H 1 1 1 1 1 1 

10 -10 1Q -9 1Q -8 1Q -7 1Q -6 1Q -5 1Q -4 10 -10 1Q -9 1Q -8 1Q -7 1Q -6 1Q -5 1Q -4 



initial ligand concentration p(M) initial ligand concentration p(M) 

Figure 3. Equilibrium time, T N) with 7V= 10, using the deterministic model of §5.2. (a) Double positive thymocytes, (b) Single 
positive thymocytes. Time units are seconds. The physiologically relevant range of initial ligand concentration is shaded grey. 



(For example, compare panels a and b (change in r) and 
panels a and d (change in N) in figure 1.) 

S.S. Equilibrium time to N TCR-pMHC 
complexes 

If one was to assume that ligand hierarchy or potency (for 
a given TCR) is determined by how rapidly thermal 
equilibrium of TCR-pMHC complexes is established, it 
would be natural to compare the dose- response binding 
curves for the different ligands under consideration [19]. 
In figure 2, we plot f eq (see §5.2), the fraction of free 
TCR molecules at equilibrium, as a function of the ini- 
tial ligand concentration (dose-response) for different 
ligands and for DPs (a) and SPs (b). These curves, 
which have been computed making use of the equations 
derived in §5.2, agree with the experimental equilibrium 
binding curves of Naeher et al. [19]. 

We now continue to explore the equilibrium 
dynamics of receptor -ligand binding with a second 
T cell response criterion. We introduce T N: the time 
needed to reach N TCR-pMHC complexes. In §5.2, 
we derive an expression for T N from the solution of an 
ordinary differential equation [44]. This criterion is, to 
some extent, the deterministic version of the stochastic 
criterion that we have introduced in the previous 



section. In figure 3, we show the time to reach N= 10 
TCR-pMHC complexes for both double and single 
positive thymocytes (figure 3a and 6, respectively). 

Inspection of both the DP and SP thymocytes cases 
(figure 3a and 6, respectively), clearly indicates that 
this second equilibrium criterion does not allow dis- 
crimination between ligands for intermediate to high 
initial ligand concentrations. That is, for a large range 
of initial concentrations, this criterion cannot dis- 
tinguish between positively and negatively selecting 
ligands, contradicting what has been experimentally 
found in Naeher et al. [19]. Thus, this criterion cannot 
account for the hierarchy of activity of peptides 
observed in vivo [19]. Furthermore, for low enough con- 
centrations, there is an abrupt change in the behaviour 
of the ligands 4N and 4A, as T N becomes unbound, that 
is, for these two ligands, T N becomes uncontrollably 
large, and for small initial concentrations these two 
ligands will not reach N TCR-pMHC complexes in a 
finite amount of time. This behaviour shifts towards 
higher concentrations as the value of N increases. The 
immunological implications of this criterion are: for 
high concentrations, all the ligands elicit the same 
response and thus T N behaves as an all-or-nothing 
T cell response mechanism. If we were to include the 
condition that each TCR-pMHC complex needs to 
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remain bound for a minimum dwell time (defined later 
as t), this would only shift the curves by an amount 
r vertically. 

In the case of DP thymocytes, figure 3a, the impli- 
cations of this equilibrium criterion are even more 
disappointing as the roles of 4P and 4A are reversed and 
the threshold ligand 4A becomes the negatively selecting 
ligand. Both implications are in disagreement with the 
hierarchy of ligands experimentally determined [19,20]. 

From an immunological perspective, and as a function 
of the initial concentration of ligand, one expects the 
following behaviour [19,20]: (i) for large enough con- 
centrations (above the physiological range), negatively 
selecting ligands remain negatively selecting, positi- 
vely selecting ligands remain positively selecting and 
threshold ligands become negatively selecting, and (ii) 
for physiological concentrations, negatively selecting 
ligands remain negatively selecting, positively selecting 
ligands remain positively selecting and threshold ligands 
remain threshold. We can conclude that this behaviour 
is well characterized by the stochastic criterion intro- 
duced in §3.2 (see, figure 16,c or figure 66,c) but is at 
odds with the equilibrium criterion discussed in this 
section. 



10 4 n 




1 ~i 1 1 1 1 1 1 1 

1Q -12 1Q -ll 10 -10 1Q -9 1Q -8 1Q -7 1Q -6 1Q -5 

initial ligand concentration p(M) 

Figure 4. Mean time to signal initiation, T(N,t), for N= 10 
and r = 4 s, as a function of the initial ligand concentration. 
The ligand is 4P and the T cells are SPs. The solid line rep- 
resents T(N,t), calculated with the analytical expression 
(3.1), and the shaded area is one standard deviation, calcu- 
lated with the analytical expression (3.2). The open circles 
represent T(N,t), numerically calculated with just one 
computational run (see §5.4). Time units are seconds. 



3.4- Statistics of the mean time to signal 
initiation 

Here, we introduce the equation for the time it takes a 
T cell to reach (for the first time and not necessarily in a 
simultaneous fashion) N TCR engagements with 
pMHC ligands, such that each engagement lasts at 
least a dwell time r. We also explore the implications 
of this equation for the MTSI (which is in the theory 
of stochastic processes a mean FPT [58]), T(7V,r), as a 
function of its relevant parameters (at a given tempera- 
ture). The relevant parameters are the initial pMHC 
ligand concentration, p, the average number of TCRs 
expressed on the surface of a T cell, 7V R , and the associ- 
ation and dissociation rates for a given pMHC ligand. 
As derived in §5.3, the MTSI is given by 



T(N,t) = r + 



Ne rk oii 



(3.1) 



The previous equation can be intuitively understood as 
follows. The probability that, once engaged, a receptor 
stays engaged for at least a time r is exp(— r£; 0 ff). The 
mean number of complexes formed such that N of 
them remain bound for longer than r is N' = exp(rA; 0 ff )N. 
Note that T(N,t) is always a decreasing function of p 
and & on , but an increasing function of k Q ^. In the limit 
of very small initial ligand concentration, we have 
T(N, r) oc 1/p. In the limit of large initial ligand con- 
centration, we have T(N : r) — > r. In §5, we have also 
derived the following analytical expression for the 
variance of the MTSI, as a function of N and r: 



var(MTSI) = 



Ne 



T& 0 f 



2 • 



(3.2) 



The predictions of (3.1) are in excellent agreement with 
the exact numerical simulations (see §5) as can be seen 
in figure 4, which provides a comparison between the 



numerical simulations and the theoretical predictions 
of (3.1) and (3.2). 

In figure 5, we explore the dependence of T(7V,t), 
given an initial ligand concentration, as a function of 
t (for fixed N) and as a function of N (for fixed r). 
The left panel ( figure 5 a) illustrates the dramatic differ- 
ence in T(N,t) based on small differences in A; off . As 
predicted by equation (3.1), T(N,t) depends linearly 
on N ( figure 5 b) . 

We conclude this section with a final application of the 
stochastic criterion to the problem of ligand self- non-self 
discrimination. In the spirit of the example discussed 
by van der Merwe & Dushek [35] (box 1), we study the 
implications of the MTSI formula to illustrate the possi- 
bility of self- non-self discrimination. Suppose that 
ligand F 'foreign' has A; Q ff = 1.0 s -1 and ligand S 'self 
has k Q ff = 5.0 s -1 , that both have the same value of A^ n , 
but the concentration of S is 100 times that of F. If 
k on N R p = 10 s" 1 for F, and if k on N R p = 1000 s" 1 for S, 
then T(10, 4) = 59 s for F and T(10, 4) = 5 x 10 6 s for 
S. In this example, the stochastic criterion produces 
very clear self- non-self discrimination. 



3.5. Pre- selection DPs versus SPs 

We now explore the MTSI hypothesis on pre-selection 
DP thymocytes, which have 10-fold lower average 
number of TCRs on their surface than SP thymocytes. 
We also note that the binding rates (for a fixed temp- 
erature) are different for SP and DP thymocytes 
(table 1). We have made use of the three different 
ligands (4P, 4A and 4N) of the Tl TCR system [19]. 
Our results are summarized in figure 6. DPs are not 
as capable of discriminating between negatively select- 
ing and threshold ligands as SPs at low initial pMHC 
ligand concentrations. 
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Figure 5. Mean time to signal initiation according to the stochastic criterion. MTSI, T(N,t), for single positive thymocytes (SP), 
Tl TCRs at 37 degrees, (a) MTSI as a function of r with N = 10 and (b) MTSI as a function of TV with r= 4 s. (a, 6) The initial 
ligand concentration has been set equal to p= 10~ 8 M. Time units are seconds. 




i 1 1 1 1 1 1 ~~\ 1 1 1 1 1 1 i i 1 1 1 1 1 1 

10 -io 10 -9 10 -s 10 -7 10 -6 10 -5 10 -4 io- 10 io- 9 io- 8 io- 7 io- 6 io- 5 io- 4 io- 10 io- 9 io- 8 io- 7 io- 6 io- 5 io- 4 



initial ligand concentration p(M) initial ligand concentration p(M) initial ligand concentration p(M) 

Figure 6. Mean time to signal initiation (MTSI) for a T cell according to the stochastic criterion. MTSI, T(iV, r), for double positive 
thymocytes (DP), Tl TCRs at 37 degrees. Different panels stand for different pairs of values (N, t): (a) for (10, 1), (b) for (10, 4), (c) 
for (10, 8), (d) for (100, 1), (e) for (100, 4) and (/) for (100, 8). Time units are seconds. The physiologically relevant range of initial 
ligand concentration is shaded grey, and the dotted lines correspond to a time scale of 1 min, 1 h and 1 day, respectively. 



As can be seen from plots (a) and (cf), with r = 1 s, 
the stochastic criterion does not distinguish between 
4P (negatively selecting) and 4A (threshold). The ability 
to discriminate at 37 degrees improves if we choose r = 
4 s, (figure 66,e). These results suggest that if negative 
selection happens too early (at the pre-selection DP 
stage), it might not be effective, as pre-selection DPs 
will not distinguish between signals delivered by posi- 
tively selecting ligands and negatively selecting ligands. 
Threshold ligands, such as 4 A, are expected to behave 
as positively selecting ligands at low concentration, but 
as negatively selecting ligands at high concentration, as 
experimentally shown [19,20]. 



3.6. Two-dimensional binding data: stochastic 
T cell response criterion 

Single-molecule microscopy and fluorescence resonance 
energy transfer [55] , as well as adhesion frequency and ther- 
mal fluctuations assays [56] are recent alternative methods 
of quantifying TCR-pMHC binding, with the advan- 
tage that both molecules are anchored on cell surfaces 
[35,55,56,60]. Recent experimental kinetic parameters 
(such as k on and & 0 ff), measured in such two-dimensional 
conditions [55,56], are given in tables 2 and 6. 

In this section, we explore the stochastic criterion in 
light of these recent two-dimensional measurements. 
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Table 2. A summary of in situ (two-dimensional) binding data taken from Huppa et al. [55] from the experimental CD4 + 
TCR mouse model. 



2B4 



Ligand 



experiment 



K d (pM) 



k/2 (s) 



K (s' 



k>& = k- (s x ) 



IE k /MCC 24°C 
IE k /MCC 28°C 

IE k /MCC 33°C 
IE k /MCC 37°C 



m situ 
in situ 
in situ 



4.26 
4.33 
4.02 
4.80 



1.68 
1.25 
0.405 
0.109 



3.2 x 10" 10 

4.3 x 10" 10 

1.4 x 10" 9 
4.4 x 10" 9 



4.1 x 10" 1 
5.5 x 10" 1 
1.7 
6.4 



Table 3. A summary of in vitro (three-dimensional) binding data taken from Huppa et al. [55] from the experimental CD4 + 
2B4 TCR mouse model. 



ligand 


experiment 


K d (jjlM) 


k/2 (s) 




kos = k- (s *) 


IE k /MCC 24°C 
IE k /MCC 28°C 

IE k /MCC 33°C 
IE k /MCC 37°C 


in vitro 
in vitro 
in vitro 
in vitro 


22.9 
24.7 
27.7 
39.8 


5.77 
4.0 
2.28 
1.24 


1.7 x 10" 11 
2.3 x 10" 11 
3.7 x 10" 11 
4.7 x 10" 11 


1.2 x 10" 1 
1.7 x 10" 1 
3.0 x 10" 1 
5.6 x 10" 1 



Note that the stochastic criterion allows us to consider 
both scenarios (two- and three-dimensional), as follows. 

(i) For two-dimensional binding (receptor and ligand 
molecules anchored on cell membranes), with 
number densities of receptor and ligand M R , M L , 
respectively, A c the area of contact between the 
cells, and A;^ n wo " dimensional) and ^ t f f°" dimensional) the 
two-dimensional on and off rates, respectively, 
the MTSI is given by 



(two-dimensional) ^jy ^ 

,(tw- " 

7Ve T %ff 



(3.3) 



^(two-dimensional) ^ cMrMl ' 



(ii) For three-dimensional binding (receptor mole- 
cules anchored on cell membrane and ligand 
in solution) , with ligand concentration p, number 
of receptors on the cell surface N R and 

^(three-dimensional) ^ ^three-dimensional) ^ three _ 

dimensional on and off rates, respectively, the 
mean MTSI is given by 

^ ( t hree-dimensional) 

^(three-dimensional) /^y- j _[_ ^ 

^ ' ' ^(three-dimensional) ^ 

(3.4) 



We make use of the stochastic criterion with exper- 
imental data, measured in a two-dimensional context, 
from earlier studies [55,56]. Let us assume that given 
a TCR, ligand potency correlates inversely with the 
value of T(N,r), that is, the most potent ligand for 
the given TCR is that with the smallest value of the 
MTSI. We note that for a chosen value of r, the 
ligand hierarchy does not depend on TV, as can be seen 
from (3.3) and (3.4). 



The experimental set-up of Huppa et al. [55] is for 
the CD4 + 2B4 TCR mouse model and for the IE k / 
MCC ligand at different temperatures. We have con- 
sidered their binding data for different temperatures 
and both in situ and in vitro conditions [55] (tables 2 
and 3). The stochastic criterion (with N= 10 and r E 
[1,10] s) has been used to rank the experimental data 
(fixed ligand at different temperatures) according to 
their MTSI values (tables 4 and 5). If we assume that 
the MTSI correlates with the time of a T cell response, 
and thus, with the potency of the ligand, 1 the stochastic 
criterion implies that it is only for r > 5 s when 
both in situ (two-dimensional) and in vitro (three- 
dimensional) binding data agree on ligand potency (or 
hierarchy). In this case, and for r>5s, it is the 
lowest temperature that yields the shortest MTSI, and 
as the temperature is increased, the MTSI increases. 
A plausible way to reconcile both in situ and in vitro 
data is to hypothesize that the early intracellular mol- 
ecular steps of the T cell signalling cascade (a kinetic 
proof-reading mechanism [52,61-63]), require at least 
a time r > 5 s. 

The experimental set-up of Huang et al. [56] is for the 
CD8 + OT1 TCR mouse model and six different ligands. 
We have considered their two-dimensional binding data 
[56] (table 6), and the three-dimensional binding data of 
Gascoigne et al. [64] (table 7). The stochastic criterion 
(with N= 10 and r E [1,10] s) has been used to rank 
the experimental data (six different ligands at a given 
temperature) according to their MTSI values (tables 8 
and 9). If we assume that the MTSI correlates with 
the time of a T cell response, and, thus, with the 
potency of the ligand, the stochastic criterion suggests 
that the two-dimensional and three-dimensional 
hierarchies are very different. The three-dimensional 
hierarchy does not change, except for a switch at r ~ 5 
s, when V-OVA and R4 interchange their rankings. On 
the other hand, the two-dimensional hierarchy depends 



1 We note that in this case, the ligand is kept fixed but the 
experimental temperature increases (tables 2 and 3). 
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Table 4. Each row is the ligand ranking, for a given value of t, according to the stochastic criterion for the parameters in 
table 2. The first ranking ligand corresponds to the shortest MTSI. The parameters and colour scheme correspond to the 
in situ experiments of table 2. 



t(s) first in ranking second in ranking third in ranking fourth in ranking 



0 IE k /MCC 37°C IE k /MCC 28°C IE k /MCC 24°C 

1 IE k /MCC 28°C IE k /MCC 24°C IE k /MCC 37°C 

2 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

3 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

4 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

5 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

6 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

7 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

8 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

9 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 

10 IE k /MCC 24°C IE k /MCC 28°C IE k /MCC 37°C 



Table 5. Each row is the ligand ranking, for a given value of t, according to the stochastic criterion for the parameters in 
table 2. The first ranking ligand corresponds to the shortest MTSI. The parameters and colour scheme correspond to the in 
vitro experiments of table 3. 



r(s) 



first in ranking 



second in ranking 



third in ranking 



fourth in ranking 



0 
1 

2 
3 
4 
5 
6 
7 
8 
9 

10 



IE k /MCC 37°C 
IE k /MCC 37°C 
IE k /MCC 37°C 



IE k /MCC 28°C 
IE k /MCC 24°C 
IE k /MCC 24°C 
IE k /MCC 24°C 
IE k /MCC 24°C 
IE k /MCC 24°C 



IE k /MCC 28°C 
IE k /MCC 28°C 
IE k /MCC 24°C 
IE k /MCC 24°C 
IE k /MCC 28°C 
IE k /MCC 28°C 
IE k /MCC 28°C 
IE k /MCC 28°C 
IE k /MCC 28°C 



IE k /MCC 28°C 
IE k /MCC 24°C 
IE k /MCC 24°C 
IE k /MCC 24°C 
IE k /MCC 28°C 



IE k /MCC 24°C 
IE k /MCC 28°C 

IE k /MCC 33°C 
IE k /MCC 37°C 
IE k /MCC 37°C 
IE k /MCC 37°C 
IE k /MCC 37°C 
IE k /MCC 37°C 
IE k /MCC 37°C 
IE k /MCC 37°C 
IE k /MCC 37°C 



Table 6. A summary of two-dimensional binding data taken from Huang et al. [56] for the CD8 + OT1 TCR experimental 
mouse model. 



ligand (25°C) 


experiment 


A c k on (^m 4 s x ) 


MO 


^off= k- (s x ) 


OVA 


two-dimensional 


1.7 x 10" 3 


7.6 x 10" 2 


7.2 


A2 


two-dimensional 


9.2 x 10" 4 


4.0 x 10" 2 


3.3 




two-dimensional 


4.7 x 10" 5 


2.1 x 10" 3 


3.4 


El 


two-dimensional 


1.1 x 10" 5 


4.9 x 10" 4 


2.6 


V-OVA 


two-dimensional 


1.7 x 10" 6 


7.6 x 10" 5 


0.9 


R4 


two-dimensional 


2.0 x 10" 6 


8.9 x 10" 5 


1.8 



a lot on the value of r. We note that the hierarchy 
of ligands for r < 1 s is almost an inversion of the hierar- 
chy that gets established for values of r greater than 6 s. 
This is not surprising, given the negative correlation 
between two-dimensional and three-dimensional off- 
rates reported by Huang et al. [56]. 



4. DISCUSSION 

T cell receptor binding to ligand determines T cell 
responses at different times and locations during the life 
of a T cell [15]. Current experimental evidence, as 
reviewed in Valitutti et al. [25], provides support to the 
following: (i) sufficiently long receptor -ligand 



engagements are required to initiate the T cell signalling 
cascade that results in productive signal transduction, 
and (ii) counting devices are at work in T cells to allow 
signal accumulation, decoding and translation into bio- 
logical responses. In other words, both the duration and 
number of TCR-pMHC bindings play a role in T cell 
responses [21,30]. These ideas have already been explored 
by Palmer & Naeher [20] in order to provide a biophysical 
basis for their affinity threshold for negative selection 
hypothesis [19]. The authors introduce the concepts of 
dwell time of individual TCR -pMHC complexes and pro- 
ductive TCR interactions, and compute the number of 
TCR-pMHC interactions required as a function of the 
TCR-pMHC complex half- life, for a given choice of 
dwell time and number of productive TCR interactions 
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Table 7. A summary of three-dimensional binding data taken from [64] for the CD8 + OT1 TCR experimental mouse model. 



ligand (25° C) 


experiment 


K d (mM) 


k/2 (s) 


k on (s^M" 1 ) 


Ms" 1 ) 




&off = k- (s x ) 


OVA 


three-dimensional 


6.2 


33 


3388 


1.1 x 10" 


n 


2.1 x 10" 2 


A2 


three-dimensional 


4.2 


34.7 


4756 


1.6 x 10" 


n 


2.0 x 10" 2 




three-dimensional 


10 


99 


700 


2.3 x 10" 


12 


7.0 x 10" 3 


El 


three-dimensional 


20.8 


10.3 


3235 


1.1 x 10" 


11 


6.7 x 10" 2 


V-OVA 


three-dimensional 


22.9 


18.7 


1619 


5.4 x 10" 


12 


3.7 x 10" 2 


R4 


three-dimensional 


48.9 


5.4 


2625 


8.8 x 10" 


12 


1.3 x 10" 1 



Table 8. Each row is the ligand ranking, for a given value of t, according to the stochastic criterion for the parameters in 
table 6. The first ranking ligand corresponds to the shortest MTSI. The parameters and colour scheme correspond to the two- 
dimensional experiments of table 6. 



t (s) first in ranking second in ranking third in ranking fourth in ranking fifth in ranking sixth in ranking 



0 


OVA 


A2 




El 


R4 


V-OVA 


1 


A2 




OVA 


El 


V-OVA 


R4 


2 


A2 


V-OVA 


El 


R4 




OVA 


3 


V-OVA 


A2 


R4 


El 




OVA 


4 


V-OVA 


A2 


R4 


El 




OVA 


5 


V-OVA 


R4 


A2 


El 




OVA 


6 


V-OVA 


R4 


A2 


El 




OVA 


7 


V-OVA 


R4 


El 


A2 




OVA 


8 


V-OVA 


R4 


El 


A2 




OVA 


9 


V-OVA 


R4 


El 


A2 




OVA 


10 


V-OVA 


R4 


El 


A2 




OVA 



Table 9 


Each row is the ligand ranking, for 


a given value of t, 


according to the stochastic criterion for the parameters in 


table 7. 


The first ranking 


ligand corresponds 


to the shortest MTSI. The parameters and colour scheme correspond to the 


three-dimensional experiments of table 7. 










t(s) 


first in ranking 


second in ranking 


third in ranking 


fourth in ranking 


fifth in ranking 


sixth in ranking 


0 


A2 


OVA 


El 


R4 


V-OVA 




1 


A2 


OVA 


El 


R4 


V-OVA 




2 


A2 


OVA 


El 


R4 


V-OVA 




3 


A2 


OVA 


El 


R4 


V-OVA 




4 


A2 


OVA 


El 


R4 


V-OVA 




5 


A2 


OVA 


El 


R4 


V-OVA 




6 


A2 


OVA 


El 


V-OVA 


R4 




7 


A2 


OVA 


El 


V-OVA 


R4 




8 


A2 


OVA 


El 


V-OVA 


R4 




9 


A2 


OVA 


El 


V-OVA 


R4 




10 


A2 


OVA 


El 


V-OVA 


R4 





(see fig. 3 of Palmer & Naeher [20]). In light of these 
results, and the fact that in the thymus SP thymocytes 
only have 4-5 days to scan the medullary environment 
[12,65,66] and in the periphery the dose- and time-depen- 
dence of antigen localization determine whether 
protective immunity is induced or not [67], in this 
study we have explored, with mathematical models, the 
timescales associated with T cell responses. We have con- 
sidered two different criteria: a stochastic one — (i) the 
mean time (or mean FPT) it takes to have had TV recep- 
tor- ligand complexes bound for at least a dwell time, r, 
each — and one based on equilibrium — (ii) the time a 
threshold number, TV, of TCRs are occupied at equili- 
brium. The dynamics of a small number of TCR- 
pMHC binding events is naturally described as a stochas- 
tic process, without the need to assume that TCR- 
pMHC association/dissociation kinetics has reached 



thermal equilibrium [36-40]. We have applied both the 
deterministic and stochastic criteria to the experimental 
data presented in Naeher et al. [19]. Our results indicate 
that the stochastic criterion provides support to the 
thymic affinity threshold hypothesis suggested by 
Palmer [19,20], whereas the equilibrium one does not. 
Thus, we propose as a timescale associated to T cell 
responses: the first time at which the stochastic criterion 
is satisfied, which is referred to as an FPT in the theory of 
stochastic processes [58] , but has been referred to as the 
first time to signal initiation in this study. Furthermore, 
other properties of the stochastic criterion are (i) for 
the values of N and r considered, the calculated MTSIs 
are of the order of the timescales associated to negative 
selection and T cell activation, (ii) a very small number 
(fewer than 10) of cognate ligand molecules is sufficient 
to elicit a T cell response, which is consistent with the 
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serial engagement model, (iii) it provides an intuitive way 
to understand self- non-self discrimination, (iv) it rele- 
gates TV to a secondary role to that of the dwell time, t, 
which is consistent with the kinetic proof-reading model, 
and (v) it can be applied to either two- or three-dimen- 
sional binding data. Our results indicate that for the 
experimental data of Huppa et al. [55] one can identify 
a threshold dwell time, r, for TCR-pMHC complexes 
of 5 s that can account for the same ligand hierarchy 
for both in situ and in vitro data. For the data reported 
by Huang et al. [56] , we note that the two-dimensional 
hierarchy, as determined by the stochastic criterion, is 
extremely sensitive to the value of r: the hierarchy of 
ligands for r < 1 s is almost an inversion of the hierarchy 
that gets established for values of r greater than 
6 s. This is not surprising, given the negative correla- 
tion between two-dimensional and three-dimensional 
off-rates reported by the authors. 

The previous discussion also stresses one of the main 
results of this study, namely that stochastic effects are 
important in the timescales that determine cellular 
responses: during thymic negative selection [12,65,66] 
or during the initiation of T cell responses in the 
periphery [67]. 

As discussed in Altan-Bonnet & Germain [52], rapid, 
sensitive and highly discriminatory TCR-induced sig- 
nals, yet exquisitely ligand specific, can be explained in 
terms of a negative feedback (phosphatase mediated), 
which suppresses signalling by weak ligands, and a posi- 
tive feedback (ERK mediated), which is induced by 
strong TCR-pMHC ligands. Furthermore, recent work 
by Dushek et al. [63] has considered the role of TCR- 
pMHC rebinding, within a kinetic proofreading model 
[61,62], as a potential mechanism for pMHC ligand dis- 
crimination. In this study, we have not attempted to 
provide a mechanistic derivation of the dwell time, t, 
introduced in the stochastic criterion. The mechanisms 
discussed in earlier studies [52,63,68] will allow us to 
explore different kinetic proofreading scenarios and feed- 
back loops, to justify the origin of r. This and TCR/ 
pMHC diffusion on cellular membranes [69] will be 
considered elsewhere. 

We conclude by summarizing the mathematical 
results of this study: we have made use of a stochastic 
model for the binding and unbinding kinetics of recep- 
tor-ligand interactions [38], that is a birth and death 
stochastic process for the number of receptor -ligand 
complexes. This model has allowed us to formulate the 
stochastic criterion and to derive analytical expressions 
for the mean value of TSI, T(N,t) as a function of N 
and r, and its variance. 



5. METHODS 

5.1. Stochastic model 

We study the dynamics of monovalent receptor binding 
to monomeric ligand with a stochastic model that is 
represented as follows: 

* I h 1*1 



An unbound receptor can bind a free ligand with rate 
and an engaged receptor can become dissociated 
from the ligand with rate Both k ± have units of 
inverse time. At the initial time, t= 0, all M R receptors 
are unbound, and there are M L free ligands. The sto- 
chastic variable represents the number of engaged 
receptors at time t. Its state space S is given by the set 
S = {0,1,2, . ..,ram(M R ,M L )}. The dynamics of the 
stochastic variable (number of engaged receptors 
at time t) can be derived from the transition probabil- 
ities that prescribe the events that can take place in a 
small time interval. There are only two types of events: 

(i) an association event that increases the number of 
engaged receptors by one unit, and 

(ii) a dissociation event that reduces the number of 
engaged receptors by one unit. 

The stochastic model for the binding and unbinding of 
receptors and ligands is a continuous time Markov 
chain, a birth and death process [58,70] with rates 

H n = k-n, 0 < n < min(MR, Ml), 

X n = k + (M R - n)(M L - n), 0<n< min(M R , M L ). 

Let p n (t) be the probability that there are n engaged 
receptors at time t. That is, p n (t) = ProbpQ = n]. 
These probabilities satisfy the following system of 
differential equations [70]: 



d 
dt 
d 
dt 



po = ia x pi - A 0 po, (5.1) 

Pn = + K-lPn-1 ~ (fl n + K)Pm 



1< n < M R - 1 



(5.2) 



and — Pm h ^m r -iPm r -i M*m r Pm r , (5.3) 

where we have assumed that M R < M L . Similar 
equations can be derived in the case M R > M L . The 
mean number of engaged receptors at time t is 
*(*) = E(X t ) = E^o n P n(t). 

It obeys the following differential equation: 

d a d 

— x= > n — p n 
dt ^ dr 

n=0 
M R -1 

= ^2 n [^n+lPn+l + K-lPn-l ~ (fl n + K)Pn] 

n=l 

+ M r (Am r -iPm r -i - ^m r Pm r ) 

Mr Mr-1 

= - Yl Pn + Yl XnPn 

n—1 n=0 



Mr 



Mr 



-*L. J2 n Pn + K ( M R " n )( M L - n )Pn 



n=0 



-k_x + k+M L M K - k+(M L + M n )x + k+y, 



(5.4) 



where y(t) = E(Xj)=^ 0 n 2 Pn(t) 
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If the experimental system under consideration has a 
number, 7V C , of T cells, each of them with an average 
number, 7V R , of TCRs on their surface, the mean number 
of engaged receptors per T cell at time t is given by 



(5.5) 



Let us introduce the stochastic variable, Z^, such that = 
N c Z t , and the variable w(t) = E(Z^). If we make use of 
these definitions and the fact that M R = N C N R , it is easy 
to derive the following equation for the time evolution of 
z(t): 



dt 



z = -k-z + k+M L N K - k+(M L + N c N K )z 
+ k+N,w. 



(5.6) 



5.2. Deterministic approximation 

The deterministic approximation consists of neglecting 
fluctuations in equation (5.6) by setting w(t) = z 2 (t), 
so that dz/dt = f(z), where 

f(z) = k + M L N R - [k- + h(M L + N c N R )]z + k + N c z 2 



-k-z + k+(N R - z)(M L - N c z). 



The two solutions of f(z) = 0 are z 1 and z 2 , where 



Z2 



k_ + k+(M L + N c N R ) 



'k_ + k + (M L + N c N R ) 



k+N c 



- 8 
+ 5 



and 



8 = 



[Al. + K(M L + N C N R )} 2 - 4k 2 + M L N c N R 



(5.7) 

(5.8) 
(5.9) 



(5.10) 



As t — > oo, z(t) — > z eq = z\, the stable steady-state 
value [44]. We also introduce the per T cell fraction of 
unbound receptors in the steady-state, given by 
/ eq = (N R — Zi)/N R . The exact solution of equation 
(5.6), z(t), with initial conditions, z(0) = 0, is given by 



z(t) z x z 2 



1 



-k+8t 



Z2 ~ Zi& 



-k+St 



An important quantity that is obtained from this sol- 
ution is the time to reach N TCR-pMHC complexes, 
that is, TV, such that z( T N ) = N. T N can be computed 
to yield: 



5.2.1. Ordinary differential equation under the assump- 
tion of soluble ligand binding 
We may rewrite equation (5.6) as 



dt 



z — koffZ + k 0 



N R p- z( p + ^tt^ ) +" 



VN A J VN A 



w 



where we have introduced the following parameters: 
p = initial concentration of ligand in the 



experimental system = 



VN A ' 



k 0 ii 



and t 



1/2 



log2 

k 0 {{ 



k 0 n = VN A k + and K d 



kon 



with V the volume of the experiment and N A 
Avogadro's number. 

The deterministic approximation consists of neglect- 
ing fluctuations by setting y(t) = x 2 (t), or w(t) = z 2 (t), 
so that 



dt 

If we introduce 

N C N K Mr 



kottZ+ko n (N R - z) [p 



VN AP M L ' A± 2 



VN A 



(5.11) 



+ < 



1+^(1-^7) 

lid 



4p 



and A = A + — A_ , 



then the solution of equation (5.6), with z(0) = 0, is 
given by 



z(t) = N R 



k + k-K d (1 - e - Mofft ) 
prj (A+ - X-e~ Ak o^) ' 



(5.12) 



From this solution, one can obtain the equilibrium 
value of the average number of engaged receptors per 
T cell, £ eq , and T N . These are given by the following 
expressions: 



and 



lim z(t) z eq N R — — . 



1 (N R k + X_K d -NX_p V 
-log' 



AbffA *\N R k+k-K d - NX+pi] 



(5.13) 



(5.14) 



5.3. The mean and variance of the first passage 
time {or first time to signal initiation) 

We now suppose that T cell responses take place once N 
TCRs have been engaged with ligand, for at least a time 
r each. The first time at which this criterion is satisfied 
is referred to as a FPT [58]. We will derive expressions 
for its mean value, T(N,r), and variance. 

At the instant of its formation, any ligand -receptor 
complex has probability exp(— rA: 0 ff) of remaining bound 
for longer than the dwell time r. A binding that does so is 
said to be productive. Let N f be the mean total number 
of binding events before the Mh productive one. We can 
then write (as binding events are independent) 

N f = exp(r£; 0 ff)7V. 

Let ti be the time that the ith ligand -receptor complex is 
formed. By definition, we have that FPT (TV, r) = r+ t N > 
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and thus, T = t+ Each of the TV times between 

binding events, t i+ i — t„ is exponentially distributed and 
therefore has standard deviation equal to its mean [58]. 
The time t N >, which is a sum of exponentially distributed 
random times, therefore has variance proportional to N' . 
If, up to time t N >, the number of bound receptors is much 
less than N R and much less than pVN A /N c1 then t N > is 
the sum of N' independent, exponentially distributed 
random variables with mean (& 0 nP^R)~\ that is 



N' 



and var(^) 



N' 



-, (5.15) 



so that 



E(FPT) = T(N,t) = t + 



and 



var(FPT) = var(r+ t N >) = var(fo/) = 



7Ve T 



(5.16) 



Note that the coefficient of variation of the FPT can be 
computed as follows: 



standard deviation (FPT) 
mean (FPT) 



1 



/W^Tk on pN R /VW 

(5.17) 



With more generality, we may assume that the initial 
concentration of ligands, p, is constant (neglect ligand 
depletion due to binding), but take into account that 
the number of occupied TCRs on a T cell is a finite frac- 
tion of N R . The total number of receptors that bind at 
least once before time t is N R [1 — exp(— k on pt)]. Let 
B(t) be the mean number of bound receptors at time 
t. Then, we can write 



— B= k on p(N R - B)-k oE B, 

and if we solve for B(t), we have 
k on p 



(5.18) 



B(t) = N R 



k>nP + k>B 



[l-eM-(konP+k oS )t)}. (5.19) 



The mean number of binding events up to time t is C(t), 
where 



p(N R -B(s))ds 




o 



B(s) ds 



t + 



k on p 



(5.20) 



k 0 nP+k oii (k on p+k 0 K) 2 



x(l -exp(-(^ n p-h koz)t))j. 

If (k on p + k oii )t < 1, then C(t) = N R k on p(t - \ 
k on pt 2 + •••)• Setting C(T — r) = N gives the following 



expression for T: 

T(N,t) = r + - 



AT 1 / AT 
^ + 2 V^Vr 



+ ■ 



which includes correction terms to the approximation of 
T(N, r) calculated in (5.16). Note that the correction 
terms become negligible when N R ^> N. For the 
values of N considered in this study, the correction 
terms can be neglected, as has been verified in numeri- 
cal computations. 

We conclude by providing, without derivation, the 
analytical expressions for the mean and the variance 
of the FPT in a two-dimensional setting. It is easy to 
show that 



(two-dimensional) 



(FPT) = T, 



(two-dimensional) 



(N,t) 



, (two-dimensional) 

TVe^oft 



(two-dimensional) 



"on 



A C M R M L 



and 



var. 



(two-dimensional) 



(FPT) = 



, (two-dimensional) 

7Ve T ^off 



(^on 



(two-dimensional) a n t J/f \2" 
A c Mr Ml ) 

(5.21) 



5.4- Numerical simulation method 

We use the stochastic simulation, or Gillespie algorithm 
[71,72], where the number of bound ligands as a function 
of time, in each realization, is explicitly generated. If 
there are n bound ligands at time t, then the first event 
after time t is either binding, with probability A n / (A n + 
fin) , or unbinding, with probability [x n / ( A n + /x n ) . The 
time at which the event occurs is t = At, where At is a 
random variable and Prob[A£ > s] = exp(— (A n + aO 5 )? 
s > 0. We study the dynamics of the receptor -ligand 
system using parameters N R and M L corresponding 

to the ligands 4P, 4 A and 4N. Binding and unbinding times 
are recorded, and a realization ends when the stochastic cri- 
terion (TV bindings that last, for at least, a time reach) is 
first satisfied. By averaging over realizations, we are able 
to compute mean FPTs and coefficients of variation. 
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